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. . Abstract 

Long-lived stau shows up in various super symmetric models, like gauge-mediated 
supersymmetry (SUSY) breaking model. At the LHC experiment, long-lived stau is 
useful not only for the discovery of SUSY signals but also for the study of the detailed 
properties of superparticles. We discuss a method to obtain information on spins and 
chiralities of superparticles in the framework with long-lived stau. We also show that 
such a study can be used to distinguish SUSY model from other models of new physics, 
like the universal extra dimension model. 



Low energy supersymmetry (SUSY) is an attractive candidate of the physics beyond the 
standard model, and superparticles are important targets of the LHC experiment. Even 
if signals of superparticles are found, however, it is non-trivial to confirm that the newly 
discovered particles are superparticles. This is because the SUSY-like mass spectrum is pos- 
sible in some class of models other than supersymmetric one. Thus, once exotic particles are 
found at the LHC, their properties should be studied in detail to understand the underlying 
model. 

Procedure to study the properties of superparticles crucially depends on their mass spec- 
trum; in particular, for each candidate of the lightest supersymmetric particle (LSP) in 
the minimal supersymmetric standard model (MSSM) sector, which we call MSSM-LSP, we 
expect different type of SUSY signals at the LHC experiment. Even though the lightest 
neutralino is the most popular candidate of the MSSM-LSP, charged (and/or colored) su- 
perparticle can also be the MSSM-LSP if it is unstable. Thus, for each candidate of the 
MSSM-LSP, we should consider how and how well the SUSY events can be studied at the 
LHC. 

In the present study, we consider an important possibility that the lighter stau f is the 
MSSM-LSP. Stau can be the MSSM-LSP in well-motivated SUSY breaking scenarios, like the 
gauge mediated SUSY breaking scenario [U [2] . Even though f is expected to be unstable 
in such a case, its decay length can be much longer than the size of the LHC detectors 
(~ 10 m). Then, f is regarded as a long-lived charged particle at the LHC experiment. For 
example, in the gauge-mediation model, f decays into r + gravitino, and the decay length 
becomes longer than 10 m if the gravitino mass is heavier than ~ 1 keV. Then, in such a 
case, we can observe a track of f at the LHC experiment and we expect very unique signals 
at the LHC. The f track should be useful not only for the discovery but also for the study 
of the properties of superparticle [31 III El El El El El [121 [III [13 [131 [III [13 [IHl [13 [IHl HH| ; in 
particular, in the previous study, we have shown that the masses of squark, sleptons, and 
neutralinos can be measured with very small errors using long-lived f [PT] . 

In this letter, we extend our previous analysis to discuss how we can study the prop- 
erties, in particular, spins and chiralities (i.e., handednesses), of the superparticles in the 
decay chain. Spin and chirality measurements for the neutralino-MSSM-LSP case have been 
discussed in many literatures [201 EH [221 E31 I2H 123 ESI EH I2H]- We will see that, if f is 
the MSSM-LSP, the study becomes easier because SUSY events can be distinguished from 
standard-model backgrounds by identifying f-track, and also because full event reconstruc- 
tion is possible due to the absence of the missing momentum. Information on the particles 
in the decay chain is extracted from the invariant-mass distribution of the decay products. 
Unfortunately, the observed invariant-mass distributions are deformed from the parton-level 
predictions, which becomes the origin of systematic uncertainties in the test of underlying 
models. We propose to analyze the ratio of the numbers of two different processes with same 
event topology, from which many of the systematic uncertainties should cancel out. We will 
see that such an analysis is useful to understand the chiralities of the superpaticles in the 
decay chain. Furthermore, we can also distinguish SUSY model from other models with 
particles with different spin, like the universal extra dimension (UED) model [21], which 
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may result in a SUSY-like mass spectrum. We note here that, in order to make our points 
clear, we consider supersymmetric standard model with long-lived f. However, the proce- 
dure should work in other class of models, like the UED model with a long-lived Kalza-Klein 
(KK) charged lepton 
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Let us start with discussing the underlying model we have in mind, i.e., supersymmetric 
model with f-MSSM-LSP. If low-energy SUSY exists, squarks and gluino are particularly 
produced at the LHC In the following discussion, we consider the case where the gluino is 
heavier than squarks. (Such a mass spectrum is naturally realized, for example, in gauge 
mediation model with f-MSSM-LSP.) In such a case, the processes pp — > uu and dd have 
significant cross section compared to other processes because of large parton densities of the 
up- and down-quarks in proton. (Other types of SUSY processes also occur, but their effects 
become subdominant in the following study by imposing relevant kinematical cuts.) Then, 
once q (where, here and hereafter, q is for u and d) is produced, it may cause the decay chain 
q — > qxl, followed by Xi ~~ >* t ±; t t (with \i being the lightest neutralino); for simplicity, we 
denote such a decay chain as q — > gr ± f T . In the following, we only use the hadronic decay 
mode of r. Then, SUSY events result in f tracks and r-jets as well as energetic jets. If the 
velocity of f is small enough, SUSY events can be distinguished from the standard-model 
events by identifying f track. The discovery of such events should be a clear indication of 
the existence of a new physics beyond the standard model. In addition, with the study of the 
endpoints or peak positions of the invariant- mass distributions, information on the masses 
of the new particles (i.e., squarks and neutralinos) will be obtained [17] . 

Once the masses of newly discovered particles are determined, the next task will be to 
precisely understand the underlying model. In particular at the early stage of the LHC 
experiment, it is non-trivial to confirm that the underlying model is MSSM. In the present 
case, the existence of particles with masses of m^ and m x o as well as a long-lived charged 
particle can be experimentally confirmed. However, if the masses are the only information 
available from the experiment, it is not clear if those particles are superparticles. In addition, 
even if the underlying model is assumed to be the MSSM, chiralities of observed q and f 
are unknown. As we have mentioned, the information on the underlying model is imprinted 
in the invariant-mass distributions of the particles from the decay; in order to confirm or 
exclude a specific underlying model, one should check the consistency between observed 
invariant-mass distribution and prediction of the postulated model. In the following, we 
discuss how and how well we can perform such an analysis. 

In the present case, invariant mass of (q, r) system contains important information. For 
the study of the SUSY model, we parameterize the relevant interaction terms as 

Ant = Xi^lPl + g q ,RP R )qq* + Xi(9t,lPl + 9r, R P R )rf* + h.c, (1) 

where g g< L, g q> R, g Tt L, and g Tt n are coupling constants. (In the following, we consider the 



* 1 In the simplest UED model, KK mode of the U(l)y is the lightest KK particle (LKP). However, such a 
mass spectrum can be easily modified by introducing the brane-localized interactions. Then, the long-lived 
KK lepton may show up when the KK mode of the graviton is the LKP while the KK mode of a lepton is 
the seconddightest KK particle. 



case that the production of right-handed squarks plays an important role and also that the 
lighter stau is right-handed. Then, for the process we will study, g q> L = g T ,L = 0.) The 
invariant-mass distributions of the decay processes are given by 

1 dT^ qT+f - (9 2 q ,L9r,L + 91r9t,r)(1 ~ X q r) + (9l L 9r,R + 9JR9r,L)x q r 



2±q^qT+f- dx qT 9q,L9 T ,L + 9q,R9 T ,R + 9q t L9 T ,R + 9q t R9 T ,L 

1 dT^ qT - f+ {gl L 9r,L + 9l R 9r, R )x q r + {9lL9r,R + 9^9^) i 1 



(3) 



2T^ qT - f+ dx qT 9q, L 9 2 r,L + 9q,R9r,R + 9q, L 9r,R + 9^9^ 

and dT q *^. qT ± f ^/dx qT = dT q ^. qT ^^±/dx qr , where 

x qT = MU Ml = (p q + VrfjMl, (4) 

with Ml being the maximal value of M qT : 

(m? - m 2 Q ){m 2 Q - m?) 
Ml = — « ^^ . 5 

Because of the Majorana nature of x?> ^q~>qr+f- = r<?->gi—f+- One can easily see that 
dT q ^, qT ± f ^/dx qT has non-trivial dependence on x qT . We also note here that the distributions 
of M qT depend on the chiralities of q and f and that the invariant-mass distributions are dif- 
ferent for q —¥ qr + f~ and q — > qr~? + . These facts are important in the following discussion. 
The distributions given in Eqs. ([2]) and ([3]) are crucial check points of the present model. 

Now, we show how the observed invariant-mass distributions behave by using the MC 
analysis. In our study, we work in the framework of gauge-mediated model. We consider the 
situation that the MSSM-LSP is lighter stau f, which is assumed to be long-lived, and that 
the processes pp — > uu and dd have large cross sections. We adopt the following parameters: 

A = 60TeV, M mess = 900 TeV, iV 5 = 3, tan/3 = 35, sign(» = +, (6) 

where A is the ratio of the F- and A-components of the SUSY breaking field, M mess is the 
messenger scale, A^ is the number of messenger multiplets in units of 5 + 5 representation of 
577(5) grand-unified group, tan/3 is the ratio of the vacuum expectation values of two Higgs 
bosons, and /x is the SUSY invariant Higgs mass. The mass spectrum of superparticles 
is calculated by using IS A JET 7.64 [30J; the result is summarized in Table [TJ The LHC 
phenomenology of this parameter point has been studied in [T7] , which has shown that the 
masses of superparticles can be determined with relatively small uncertainties. In particular, 
m qR and m x o are measured with the accuracies of ~ 10 GeV and ~ 1 GeV, respectively, 
with the luminosity of C = 100 fb~ . In addition, m? can be also determined by combining 
time-of-flight and momentum information; the expected accuracy is ~ 0.1 GeV [8J. In our 
study, we assume that the masses of these superparticles can be well determined before the 
study of the invariant-mass distributions. 
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Mass (GeV) 
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1309.39 


u L 


1231.70 


ur 


1183.97 


d L 


1234.28 


d R 


1180.19 


h 


1082.85 


h 


1195.08 


h 


1145.24 


h 


1185.83 


V\ 


388.05 


h 


396.19 


r-i 


402.57 


v T 


383.80 


Ir 


194.39 


T\ 


148.83 


X°i 


239.52 


xl 


425.92 


xl 


508.41 


xd 


548.67 


xt 


425.45 


xf 


548.43 


h 


115.01 



Table 1: Masses of the superparticles and the lightest Higgs boson h in units of GeV. The input 
parameters are A = 60 TeV, M mcss = 900 TeV, N$ = 3, tan/3 = 35, sign(//) = +. (We use the 
top-quark mass of 171.3 GeV.) 



We have generated SUSY events for yfs = 14 TeV with HERWIG 6.510 package 
(The total cross section of the SUSY events is 669.6 fb.) In order to simulate detector effects, 
we use the PGS4 detector simulator [33J with slight modification to treat stable stau; the 
momentum resolution of f is assumed to be the same as those of muons. Following [4J , we 
assume that f with 0.4 < (3 f < 0.91 can be detected with the efficiency of 100 % with no 
standard-model background. Staus with (3f > 0.91 are assumed to be identified as muons. 

In order to use the events with the decay chain q — > qx®, followed by X\ ~~ ^ t ± t 1=f , the 
following selection cuts are applied: 

(a) At least one f with the velocity 0.4 < /3 f < 0.91. 

(b) At least one r-tagged jet j T with p T > 20 GeV. 

(c) Exactly two jets j with p? > 30 GeV. 

The requirement (c) is to eliminate the gluino production events. If there exists only one f 
with 0.4 < f3f < 0.91, the highest px muon-like object is regarded as second f because two 
staus are expected in SUSY events. Then, for all the possible combinations of (j,j T ,f), we 
perform the following study. We first reconstruct the tau momentum p T assuming that tau 
and stau are from the decay of xl (whose mass is expected to be already known). Because the 
tau from the neutralino decay is highly boosted, we approximate that the three-momentum 
of r is parallel to that of j T . Then, we obtain 

PT = Zq 1 Pjr, ( 7 ) 

where 

H = "^ (8) 

m o — mi 
x? T 

Combinations with Zq > 1 is eliminated. Then, we calculate 



M q =^(p J+ p f + z^p jT y. (9) 

Using the fact that there exists a very sharp peak in the distribution of Mq, which is from 
qR production, only the combinations with m^ pca — 40 GeV < Mq < ?7u pca + 40 GeV are 



adopted, where rn?^ 1 = 1170 GeV is the position of the peak [T?]j^ 2 l We calculate the 



distribution of the following variable 



where 



x jr = M%IM] T = (p 3 + z^ PJT Y/M* T , (10) 



(Ml-m 2 )(m 2 -m 2 f ) 

M 2 T = — Xl 9 Xl . (11) 
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Figure 1: Distributions of Xj T in the reconstructed decay processes, qn — > qT + f^ (left) and 
q~n — > qr~fft (right), where Xj T is a squared invariant mass normalized by its maximal value. Here 
we take C = 100 fb _ . The theoretical predictions are also shown in the SUSY case (red dashed 
line), and for similar processes in the UED case (blue dotted line). 



The charges of j T and f, which are both observable, should be opposite for signal events. 

In Fig. [U we show the distributions of Xj T = M^ T jM^ T for (j,j+,f~) and (j,j~,f + ) 
events. In the same figure, we also plot the theoretical differential decay rates given in Eqs. 
02]) and (EJ); the normalization is determined so that the total number of events agrees with 
the MC data. The MC results and the theoretical predictions seem to agree at qualitative 
level. 

At quantitative level, however, the agreement is not perfect. We can see that the numbers 
of events are suppressed when Xj T — > and 1. These can be understood by considering the 
event configurations in those limits. When Xj T — > 0, r is emitted in the opposite direction 
to Xi m the rest frame of q. Because most of the squarks are not significantly boosted 
with the present choice of the squark masses, r in events with small Xj T are likely to have 
small pt- Because the px of j T is required to be larger than 20 GeV, the acceptance of the 
signal is suppressed for Xj T ~ 0. On the contrary, when Xj T — > 1, the momenta of r and f 
become parallel in the rest frame of q. Such events are eliminated by the isolation cut for 
the r-jet; in the present analysis, no extra activity is allowed in the cone (with the size of 
AR = 0.5) around r-jet. To see the validity of these arguments, we estimate the efficiencies 
of corresponding kinematical cuts using the momentum information on the decay products 
obtained from HERWIG output. The results are shown in Tables [2] and [31 We can see that 
the efficiencies behave as we have expected. It is also notable that the efficiencies depend on 
x 



JT- 



but are insensitive to the charges of final-state r and f. 



Another important reason of the disagreement should be the contamination of back- 
ground, which can be from fake r-jets (mis-identified QCD jets) as well as from wrong 



# 2 The peak position in the present study is found to be smaller than the input value of the squark mass by 
~ 10 GeV. This is expected to be due to the energy leakage in the jet reconstruction, and may be corrected 
once the jet energy is well calibrated. 



Bin 


UJ?,t ) 


U,3t,t + ) 


< x qr < 0.1 


0.09 


0.13 


0.1 <x qr < 0.2 


0.19 


0.21 


0.2 < x qr < 0.3 


0.24 


0.24 


0.3 < x qr < 0.4 


0.25 


0.25 


0.4 < x qT < 0.5 


0.25 


0.26 


0.5 < x qT < 0.6 


0.27 


0.27 


0.6 < x qT < 0.7 


0.27 


0.26 


0.7 < x qT < 0.8 


0.29 


0.28 


0.8 < x qT < 0.9 


0.28 


0.28 


0.9 < x qT < 1.0 


0.32 


0.28 



Table 2: Ratio of the total number of squark-decay events and that with pt(Jt) > 20 GeV. 



Bin 


UJ?>t ) 


(j,Jr,f + ) 


< x qT < 0.1 


0.97 


0.98 


0.1 < x qr < 0.2 


0.95 


0.94 


0.2 < x qr < 0.3 


0.88 


0.87 


0.3 < x qr < 0.4 


0.80 


0.79 


0.4 < x qr < 0.5 


0.72 


0.71 


0.5 < x qT < 0.6 


0.62 


0.61 


0.6 < x qT < 0.7 


0.49 


0.50 


0.7 < x qT < 0.8 


0.41 


0.39 


0.8 < x qr < 0.9 


0.29 


0.29 


0.9 < x qT < 1.0 


0.26 


0.18 



Table 3: Ratio of the total number of squark-decay event and that with ARj T f > 0.5. 

combination where r and f have different parents. Once the real data will become available, 
an accurate determination of the number of backgrounds may be possible using, for example, 
events off from the squark-mass peak (i.e., the sideband). In the present MC analysis, we 
have estimated the shape of the background from the sideband samples, and the number of 
backgrounds in each bin is inferred to be approximately universal as far as Xj T is not close 
to or 1. In the following analysis, we adopt constant background in each bin. Since the 
accurate estimation of the total number of backgrounds is difficult from the sideband sam- 
ples in the present analysis, the normalization of background is treated as a free parameter 
and is determined so that the x 2 variable defined below is minimized. 

If the effects of the deformation will be well understood in future by, for example, a 
reliable MC analysis, the invariant-mass distributions may be directly used to discriminate 
underlying models. However, it is desirable to find quantities which are insensitive to the 
effects of deformation. For this purpose, we consider the ratio of the numbers of (j,j^, f~) 



and (j,j T ,t + ) events; we define 

Ll - ln NAjJ7^)' (12) 

where Ni(j,j^, f T ) are the numbers of events in i-th bin with charges of (j T , f) being (±, =p). 
The error of Lj is given by 

bL] = N^j^f-) + N^j-^y (13) 

To discuss how well the theoretical prediction is expected to agree with experimental result, 
we calculate 



it rw^ 

L '-E ( S/? ■ < 14 > 



where L\ is the theoretical prediction which is given by 

L (th) = ln Nt^(j,j^f-) + N^ 

and the summation is over the bins in the range of 0.1 < Xj T < 0.9 (with the width of 
Axj t = 0.1); in order to minimize the effects of background contamination, we do not use 
the bins at Xj T ~ and 1. Here, N^ na (j,j^, f T ) are theoretical predictions of the number 
of signal events in i-th bin calculated from Eqs. (j2]) and ([3]), while iV^ 50 -* is the number of 
background events in each bin, which is independent of i. 

Li from the MC analysis is shown in Fig. [2] for C = 100 fb _1 . In the same figure, 
theoretical predictions L\ are also shown. We can see that the MC result well agrees with 
theoretical prediction of the present SUSY model. In addition, in Table HI we show the result 
for an ideal case where the luminosity is sufficiently large. (Here, we generate the event for 
C = 20 ab _ .) The results in the table indicate that the effects of the deformation of the 
invariant-mass distribution almost cancel out by taking the ratio. 

The value of x 2 varies as we take different sets of event samples. In order to estimate the 
typical value of x 2 , we calculate \ 2 for 20 sets of MC samples (for a fixed value of luminosity 
£), and obtain averaged value of x 2 . For C = 30 and 100 fb _1 , the averaged values are 
found to be (x 2 ) = 7.6 and 8.0, respectively, which indicates a good agreement between the 
theoretical prediction and observation. If we flip the chirality of one of q or f, then the 
value of x 2 significantly increases; we obtain (x 2 ) = 15.8 and 31.5 for C = 30 and 100 ftT 1 , 
respectively. This provides important information on the particles in the decay chain; the 
result indicates that the squarks in the observed peak and the lighter stau have the same 
chirality rather than different ones. Information on the chirality of f may be obtained from 
other observables; one of the examples is the tau polarization [19]. Then, assuming SUSY 
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Figure 2: Distribution of L^. MC result for C = 100 fb _1 is cross with error bar. Theoretical 
predictions of Li obtained after minimizing \ 2 are also shown for the SUSY case (red circles) , and 
the UED case (blue triangles). 



model as the underlying model, we obtain information on the chirality of the dominantly 
produced squarks. 

Once the data from the SUSY events are collected in the LHC experiment, we can also 
exclude class of models other than SUSY model by studying the distribution of the invariant 
mass of (q, r) system. For a model-independent analysis, we parameterize the distribution 
of x qT for the decay chain resulting in r ± (and opposite-charge long-lived particle) as 



1 



dT qT ± 



6 



J- yj-zh CitX', 



<-lT 



c 2 + 6 



[1 ± d (2x qT - 1) + c 2 x qT (l - x qT )} 



(16) 



where c\ and c 2 are constants (with \c\\ < 1 and c 2 > 0). In the present SUSY model, 



c 



(SUSY) (SUSY) 
Co 



(—1,0). If the underlying model is assumed to be the UED model, we 
would identify the particles with the masses of rriq a and m x o with KK modes of q (denoted 

as g( KK )) and neutral gauge boson (denoted as B^ '), respectively, as well as the long-lived 
charged particle with the KK mode of r (denoted as t 1 - 1 ^). In this framework, there exists 
a similar decay chain as in the present SUSY model (i.e., g( KK ) — > qB^ , followed by 



Bin U 5Lf m) Lf h ' tiUtiY} gg^T 

< x JT < 0.1 0.97 0.45 0.91 0.26 

0.1 < x JT < 0.2 0.72 0.30 0.69 0.15 

0.2 < x JT < 0.3 0.50 0.27 0.48 0.09 

0.3 < x JT < 0.4 0.29 0.25 0.29 0.05 

0.4 < x JT < 0.5 0.07 0.24 0.09 0.02 

0.5 < x JT < 0.6 -0.12 0.24 -0.09 -0.02 

0.6 < x JT < 0.7 -0.31 0.24 -0.29 -0.05 

0.7 < x JT < 0.8 -0.47 0.25 -0.48 -0.09 

0.8 < x jr < 0.9 -0.65 0.26 -0.69 -0.15 

0.9 < x jr < 1.0 -0.86 0.32 -0.91 -0.26 

Table 4: Li, 5L\ (which is the error for the case of L = 100 fb _1 ), and L\ for the case of large 
enough luminosity. The error of Li for the luminosity C is given by 5Li = 5L\ /V ' C\qq with C\qq 
being the luminosity in units of 100 fb _ . 



BJ? K) -+ tt^). Then, we obtain 



,(TJED) 



2m 4 

U) 



(UED) _ _ (UED) (UED) XJ M 7 \ 



mi mi + 2m 



mi mi + 2m 4 
ir T Xi 



(18) 



where el = +1 when the fermion / is right-handed (i.e., singlet under SU(2)l) while 

€* = — 1 when / is left-handed. One can see that the distribution is different from the 
SUSY case; this is due to the fact that the distribution depends on the spins of particles in the 
decay chain. Using the mass spectrum of the present model, (q , q> ) ~ (—0.17, 5.0) or 
(0.17, 5.0), depending on the relative chirality of the KK modes of q and r. (For comparison, 
we also show the results for the UED case in Figs. [T] and [21 and Table @] for the case 
(c\ ,c 2 ) — (—0.17,5.0), which results in a smaller value of x 2 .) Furthermore, in a 
model where the invariant-mass distribution is flat in the phase space, (c[ a , c 2 a ) = (0, 0). 
In Fig. [31 we show the contours of (x 2 ) = 14.1, which correspond to 95 % CL. bounds for 
7 degrees of freedom. In the same figure, we also show the points corresponding to various 
models. We can see that the analysis based on the ratio of the numbers of (j,j^,f~) and 
(j,j~,f + ) events is useful for the test of underlying models. We note here that, on the 
contours given in Fig. [31 iVec is estimated to be smaller than the number of signal events 
in average. Thus, in the present case, constraint on c\ vs. c<i plane can be obtained without 
knowing the maximal possible number of backgrounds in detail; this is because Li is strongly 
dependent on Xj T , as shown in Table HI in the sample point used in our analysis. In the case 
with an underlying model giving rise to a weaker Xj T -dependence of L i; careful estimation of 
the maximal number of backgrounds may be necessary to discriminate underlying models. 
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Figure 3: 95 % C.L. bounds on the c\ vs. C2 plane for integrated luminosities indicated in the 
figure. The region right to the contour is excluded for each luminosity. 



In this letter, we have considered a procedure to study the properties of superparticles in 
the case where the lighter stau is long-lived. We have shown that properties of the particles 
in the decay chain are extracted from invariant-mass distributions of the decay products 
of q. In particular, we emphasize that the effects of the deformation of the invariant-mass 
distributions are largely reduced by taking the ratio of two different decay processes with 
the same event topology, which are, in the present case, q — > qxl, followed by Xi ~^ t + t~ 
and by Xi ~^ T~f + . We also note here that, in some of the model beyond the standard 
model, a number of particles (like neutralinos) decay into two different final states which 
are charge-conjugated to each other, and that the reduction of the effects of deformation by 
taking the ratio may be useful in various cases even if there is no long-lived heavy charged 
particle. 

In the present study, we have assumed that the squark(s) responsible for the peak in the 
Mq distribution has unique chirality. This should be also experimentally confirmed. One of 
the circumstantial evidences of this may be negative observation of the decay processes of 
squarks into Wino-like chargino and neutralino (because the squarks are right-handed). In 
addition, it may be also possible to reconstruct q~i production event to determine the mass 
of the left-handed squarks, from which the observed peak in the Mq distribution may be 
understood to be from q^ production. These will be discussed elsewhere [31]. 
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